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Abstract. We propose a novel approach to intrinsic decoherence without adding new 
assumptions to standard Quantum Mechanics. We generalize the Liouville equation just by requiring the 
dynamical semigroup property of time evolution and dropping the unitarity requirement. With no 
approximations and specific statistical assumptions we find a generalized Liouville equation which 
depends on two characteristic time x t and x 2 and reduces to the usual equation in the limit x t = x 2 — > 0. 
However, for Xj and x 2 arbitrarily small but finite, our equation can be written as a finite difference 
equation which predicts state reduction to the diagonal form in the energy representation. The rate of 
decoherence becomes faster at the macroscopic limit as the energy scale of the system increases. In our 
approach the evolution time appears, a posteriori, as a statistical variable as if time evolution would take 
place randomly at average intervals x 2 , each evolution having a time width Xj. A generalized Tarn 
Mandelstam inequality is derived. The relation with previous work by Milburn is discussed. The 
agreement with recent experiments on damped Rabi oscillations is described. 



INTRODUCTION 

The existence of coherent superposition of states is the basic reason of many 
paradoxical aspects of quantum mechanics. The evolution from a coherent superposition 
state to a statistical mixture is called decoherence. This is a central problem for 
measurement theory and for the classical limit of quantum mechanics at the macroscopic 
level. The Schroedinger cat, which can be in a superposition of states dead or alive, or a 
macroscopic particle which can be in a superposition of "here" and "there" are typical 
example of paradoxes whose interpretation is still controversial. A superposition state 
gives non zero off-diagonal elements of the density operator, which originate quantum 
interference and non-classical correlation effects, von Neuman postulated the reduction 
to the diagonal form as a result of a measurement. However, this reduction remains 



mysterious since it cannot be described by a unitary Hamiltonian evolution and does not 
clarify how and when state reduction take place and what is the underlying dynamical 
process. Essentially two approaches to decoherence have been proposed: the most widely 
used (1,2) takes quantum mechanics as it is appealing to dissipation due to the interaction 
of the system with the environment or equivalently with the measuring apparatus, tracing 
at the end, over the many degree of freedom of the environment. The decoherence of the 
system is usually described by Master Equations (ME) which are derived using reasonable 
statistical assumptions and specific but arbitrary models to simulate the environment and 
its interaction with the system. From a conceptual view point, it appears rather peculiar 
that one has to invoke the environment, dissipation or the measuring apparatus to 
conclude that the cat must be dead or alive and a big particle must be here or there and 
not in a superposition of these possibilities with interference between them. It is like to 
say that the particle is going to be localized here or there because there is friction, and 
that the macroscopic limit consist just in the fact that the damping becomes stronger in 
this limit. 

A minority of scientists look for a modification of standard quantum mechanics so that 
decoherence becomes intrinsic, i.e., not related to any specific model of interaction or 
entanglement with the universe. The most widely known solution has been proposed in 
ref.(3). The authors modify phenomenologically the Schroedinger equation adding 
nonlinear terms to the Hamiltonian so that the system, at Poisson distributed times, 
undergoes a sudden localization. Their model contains two unspecified parameters: the 
frequency and the spatial extent of localization. Other similar models have been proposed 
(4,5). However, a common feature of these models is that the energy of the system is not 
preserved. This looks to us a very peculiar aspect for an intrinsic theory of decoherence, 
which in our opinion should act on the coherence and not on the energy of the system. 

A scheme for intrinsic decoherence has been proposed by the author in 1983 (6) 
assuming a finite difference Liouville equation with time step x (cronon) which has been 
connected to the time energy uncertainty relation. In ref. (7) it has been shown that this 
finite difference equation is equivalent to a semigroup master equation of the Lindblad 
form. 

More recently Milburn (8) has proposed a modification of the Liouville equation 
assuming that: "the system does not evolve continuously under unitary time evolution, 
but rather on a stochastic sequence of identical unitary transformations" according to a 
Poisson distribution. However, in Milburn theory time is a parameter, as in standard 
quantum mechanics and the Poisson distribution is assumed. 

In this paper we generalize the Liouville equation without any specific statistical 
assumption and without using any specific model for the interaction with the 
environment. We require the semigroup property of time evolution, dropping the unitarity 
requirement, we obtain an expression for the density operator which, a posteriori, can be 
interpreted as if the evolution time is not a fixed parameter but a statistical variable 
whose distribution is that of the waiting time of "line theory", i.e., a T-Poisson 
distribution. Our time distribution function contains two characteristic time, Xi and X2, 
which appear naturally in the theory as scaling times. They look the analog in time of 
GRW localization in space. Precisely time evolution appears on as a random process in 
which X2 is the average time step between two evolution and Xi is the time width of each 
evolution. In terms of waiting time statistics, it is like to say that % 2 is the average interval 
between the arrival of two "clients" at the teller in a bank and Xi is the time each "client" 
spends at the teller. Therefore,Xi and x 2 have very different physical meaning and in 



general x 2 > Xi. Our time evolution law, for Xi = X2 — > gives back the Liouville equation, 
whereas for Xi and X2 arbitrarily small but finite give irreversible state reduction to the 
diagonal form in the energy representation. Assuming Xi = X2 = x sufficiently small, we 
obtain the Milburn ME (8). Unlike in ref. (3) in our treatment energy is a constant of 
motion. Our evolution equation can be written in the form of a finite difference equation 
with time step X2. Therefore , according to our formalism, one cannot give a continuos 
and instantaneous description of time evolution but only at time steps given by the cronon 
x 2 . From the finite difference equation we obtain a generalized Tarn Mandelstam 
inequality which connects Xi and X2 with the time-energy uncertainty relation. We apply 
our formalism to many examples which can be experimentally tested. In particular we 
derive an extra diffusion term in the position spread of a free particle, decoherence of a 
free particle prepared in a Schroedinger cat-like state of two different positions; we 
demonstrate the existence of an intrinsic linewidth for a single mode e.m. field. Finally, 
we describe decoherence for a spin superposition state in a magnetic field, cancellation of 
EPR correlation under proper conditions, intrinsic damping of Rabi oscillation in a two 
level system, in agreement with recent experimental results (9, 10), is predicted. 

THE FINITE DIFFERENCE LIOUVILLE EQUATION 

The unitary time evolution of a quantum system is generally described by the Liouville- 
von Neumann equation, dp/dt = -iLp(t) , where p is the density operator, L is the 

Liouvillian super-operator, Lp = (l/ft)[H,pJ, and H is the Hamiltonian. Its formal 

_.. — Ht -Ht 

solution can be written as p(t) = e p = e h p e s . Here the time t appears as a 
parameter, not an "observable" or a statistical variable as, for example, are position, 
momentum and H. In the energy basis, Hln>=Eln>, and one has: 

p = -ico p and p (t) = p (0)e" ifl>nn,t (1) 
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where co nm = {E n -E m )/h . The degenerate case can be included in this notation by 
assuming that the same states ln> belong to the same eigenstate E n . Let us define a 
generalized density operator p defined as a coarse grain time average of p(t): 

p(t) = }dt'P(t,f)p(f); t,t'>0 (2) 



where P(t,t') is a function still to be determined. In particular if P(t,t') = 8(t-t'), 
p(t) = p(t) , so that p(t) is just a generalized form of p(t) , because P is unspecified. 
Equation (2) can be written as 

p(t) = V(L,t)p(0) ; V(L,t) = }dt'P(t,t')e- lLt ' . (3) 


We now determine P(t, t') imposing the following conditions: 
i) p(t) = p + (t)>0; Trp(t) = l and ii) V(t + t') = V(t)V(t') . 

Condition (i) identifies p(t) as a density operator. Condition (ii) is the so called 
semigroup property which ensures translational invariance of the initial condition i.e., 
p(t + t') = V(t + t')p(0) = V(t)p(t') = V(t')p(t). Note that (i) and (ii) are satisfied by 



the usual Liouville operator V(t) = e lLt . We just drop the request of Unitarity 
V + = V" 1 . Condition (i), using Eq.(2) and taking the trace on both sides, leads to: 

}dt'P(t,f) = l; P(t,t')>0 (4) 



where we used Trp(t') = l. Equation (4) defines P(t,t') as a probability distribution 
function; P(t,t')dt' can be read as the probability that the random variable t takes a value 
between t' and t'+dt'. Therefore, p(t)and V(t) appear as an average value respectively 
of p(t) and V(t) . The semigroup property (ii) can be generally satisfied assuming 
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where x 2 is a scaling time. Note that also the Liouville operator V(t) = e~ lLt can be 
written in the form of Eq.(5) taking Vj = e lLt2 (unitary). We now determine the most 
general form of V, (L) compatible with the previous requirements. Let us use the T 
function integral identity, 

(A + iB)- k =jdX^— e" A V m ; (6) 

o l W 

which is valid for A > and k > 0. Let us identify the following: 

V=A + iB; k = t/x 2 ; X = t'/x, 

where Xi is a scaling time generally different from x 2 . In this way by imposing the 
consistency of Eq.(6) with Eq.(3) for all L and the normalization condition of Eq. (4), we 
obtain B = Lx ; and A = 1 . Therefore, from Eq.(3) and Eq.(5), we have: 



(l + iLxJ 
P(t, t') is the r distribution function: 
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— I^-p(0) ; p(t) = -— Hl + iLxJpCt). 



p(t) = V(t)p(0) = 7 r^-p(0) ; p(t) = -— Hl + iLxJp(t). (7) 
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(t,f>0) (8) 
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This argument can be supported by considering P(t, t')=0 for t'<0 and using the 
uniqueness of the Fourier transform. Equation (7) provides a generalized form of the 
density operator and of the time evolution law which has been obtained by requiring the 
density operator properties (i), the semigroup property (ii) and dropping the unitary 
requirement. Therefore, we are not adding new elements to the formalism of quantum 
mechanics but, on the contrary, we are reducing the basic assumptions. Taking 
Xj = x 2 — > in Eq.(7) one obtains the Liouville limit and P(t,t') = 8(t - 1') . However, as 

will be shown later, for Xi and % 2 arbitrarily small but finite, one irreversibly approaches 

the diagonal form in the energy representation. 

Equation (7) in the Liouville limit x i = x 2 — >0 recovers again the Liouville equation. 

When Xi and % 2 are finite, the second order expansion of Eq.(7) gives 

^-^p-^p (9) 

where L 2 p = lH,|H,pN. Equation (9), with Xi = % 2 , gives the well known "phase- 
destroying" ME, deduced by many authors using a reservoir interaction model (2, 11, 12) 
or specific statistical assumptions (8), and used in Quantum Non Demolition (QND) 



measurement theory (11). It is straightforward to show from Eq.(7) that p(t) obeys the 
following finite difference equation 

p(t)-p(t-x 2 ) = _. L _ (t) = _1[ H p (t) ] (10) 

Equation (10), for Xj=x 2 — >0, again gives the Liouville equation, whereas for 
Xj=x 2 =x, it reduces to an equation proposed a long time ago (6) to describe 
irreversible state reduction to the diagonal form. The difference here is that now it has 
been derived under very general assumptions and two characteristic times appear, with 
very different physical meanings. In ref. (7), it has been shown that for Xj =x 2 =x, 
Eq.(10) is equivalent to a semigroup ME of the Lindblad form. The same considerations 
apply here to Eq.(10), taking x 2 = x and substituting H with H = HXj /x 2 . In this way 
Eq.(10) is formally identical to that considered in ref. (6). 

Equation (10) shows a very important feature. If p(t)is a solution of Eq.(10), then 
f(t)p(t) is also a solution, provided f(t + x 2 ) = f(t) . Therefore, p(t)is uniquely 
determined only within the "cronon" x 2 , i-e- for time intervals t = kx 2 , with k integer. 

This fact implies a redefinition of p(0) which, in the standard description, is the density 
operator determined at some instant t=0. This is clearly an artifact, because the possibility 
of an instantaneous measurement of a complete set of observables to determine p(0) at 
the instant t = appears as a mathematical abstraction. Our finite interval description of 
time evolution, which follows from Eq.(10), appears much more realistic because p(0) 
can be interpreted as the density operator determined in a cronon interval X2. The 
evolution of Eq.(10) on the "time grid" gives p at later time intervals, k=t/X2 > 1, and can 
be parametrized as p(k)- p(k-l) = iLXjP(k) , with k>l, or equivalently 

p(k + l) = (l + iLXj) p(k) ~e~ lLT, p(k), for Xi small enough. It is as if time evolution 

occurs discontinuously in "quantum jumps" (12) spaced by intervals x 2 with each 
evolution taking place within a width Xi. Each "jump" can be approximated by a unitary 
time evolution only if Xi is small enough. Accordingly, Eq.(3) and Eq.(7), for 
t / x 2 = k > 1 , with T(k) = (k - 1) ! , can be written as: 

1 e" t7Tl (t'/x ) k ~' 
V(k) = Jdt'P(k,t')e- lLt ' ; P(k,f) = - (k _ l} ', , (ID 

where P(k, t') can be interpreted in two ways: either i) as the well known Poisson 
distribution in k or ii) as the r-distribution function in the continuous variable t'. Unlike 
previous treatments (8) we adopt the latter. Equation (11) can be interpreted in terms of 
the waiting time statistics for k independent events. According to Eq.(ll), time evolution 
is made up of random "events", i.e., unitary time "evolution". The probability density for 
k=t/X2 events to take place by a time t' is given by Eq.(l 1). In particular, X2 is the average 
interval between two "events" ( x 2 ' is the rate) and Xi is the time width of each event. 
Therefore , the "effective evolution time is given by t = (t/x 2 )x 1 . We propose an 
interpretation of our formalism in terms of a continuos measurement theory: X2" 1 is the 
observation rate, Xi is the time width of each observation. As we shall see, this 
interpretation has a precise meaning in laser- micromaser theory of ref. (13). In our 
formalism the interaction with the environment or with the measurement apparatus is 
described by two characteristic time Xi and X2 : this characterization is intrinsic and model 



independent in the sense that it does not depend on the way the measurement is carried 
out or on the detail on the measurement apparatus. Because the Hamiltonian is a constant 
of motion, our formalism can be applied, though not exclusively, to a QND measurement. 
According to this interpretation one obtains a dynamical description of von Neumann 
state reduction. In fact, in general, in the energy representation Eq.(7) becomes 

P„, m (t) = 7 ^— ^-P„. n (0) = rF J^_ _ = e -w e -^t pnm(0) (12) 

where y nm = (l/2x 2 )ln(l + (0^ m xf ) and V nm =(l/x 2 )arctg(O nm x 1 , which describes an 
irreversible evolution. Note that irreversibility is obtained also in the limit x l — > and 
x 2 — > , provided that x\ I x 2 is finite. In general, the decoherence rates y nm increase as 
%2 becomes smaller or Xi becomes larger. Assuming, for simplicity, non degeneracy, for 
E n =E m , therefore, C0h, m =O, and p nn = p nn (0) , so that the energy is a constant of motion, 

whereas for n ■£ m, p nm — > with a rate Y n , m . Therefore, p(t) — > 2^p n , n (0)|n)(n| i.e., p 

n 

approaches the stationary diagonal form. Then a pure state remains a pure state if and 
only if it is a stationary state; otherwise the system will evolve to a statistical mixture. 
Note that, at the macroscopic limit, when the energy scale becomes very large, the 
decoherence time 1/y becomes very short. If 0) nm x 1 «1 one has y nm ~co^ m xf /2x 2 and 

v nm "'^nm/''^- This is the result one obtains by the phase diffusing ME, Eq.(9). 

Therefore, this ME is valid only if CO nm x 1 « 1, for all n and m. This assumption can be 

made only if the spectrum is bounded so that it cannot be applied to a simple harmonic 
oscillator, as is usually assumed (10). We emphasize that the basic point of our treatment 
is that time appears as a statistical variable with a distribution function P(t,t') , given by 

Eq.(8), which for k=t/X2 = 1, is a simple exponential. However, for t/X2»l, P(t,t')is a 
strongly peaked function with a mean value, <t'>, and dispersion, o, given by 

<t'>=t = (t/x 2 )x 1 ; a = V<t' 2 >-<t'> 2 =x 1A /t/x 2 . (13) 

Note that <t'> does not coincide with t but with the effective evolution time t . The 
dispersion a scales as Vt , like in a diffusion process. According to the previous 
interpretation, the dispersion of o appears as the dispersion due to k=t/X2 statistically 
independent "events" times the width of each event Xi. For k=l, o assumes its minimum 
value Xi. Therefore, Xi appears as an "inner time", i.e., the intrinsic minimum uncertainty 

/ \-l/2 

of the evolution time. The relative dispersion, a/< t'>= (t/x 2 J , goes to zero as the 
number of evolution steps, t/X2, goes to infinity. Furthermore, for t/X2 » 1, P(t,t') can be 
approximated by a Gaussian in t' with mean value and dispersion given by Eq.(13). 



COMPARISON WITH PREVIOUS WORK 

The difference between our treatment and previous descriptions of intrinsic 
decoherence (8) can be summarized as follows: In ref. (8) it is assumed that i) the system 
evolves under a random sequences of identical unitary transformations, ii) the probability 
of n transformations in a time t is given by a Poisson distribution so that (using our 



notation) p(t) = £p(n,t)e~ mLT| p(0) where p(n,t) = ((t / x 2 )" / nl)e" t/T2 . In ref. (8) Xi=x 2 = 

n-0 

Y" 1 is interpreted as a "fundamental time of the universe". Summing up the above series 
one obtains p(t) = exp((e" lLT ' - l)t / x 2 )p(0) so that 

^(t) = — (e" lLT| - l)p(t) = — (e" lHT ' /s p(t)e lHT ' /s - p(t)) . (14) 



^2 ^2 



The same philosophy has been adopted in laser-micromaser theory (13) substituting e~ lLT| 
with the non unitary operator M (13). In this case, Xi is the interaction time, x 2 A is the 
atomic injection rate and Eq. (14) becomes the Scully- Lamb ME. Equation (14), with Xi 
= x 2 = Y" 1 , is the basic result of ref. (8). If one applies Eq.(14) with this assumption to 
calculate the average value of the annihilation operator of an harmonic oscillator, one 

obtains easily < a >= y(e~ ia,y - 1) < a > . From this equation Milburn (8) infers that there 

are particular frequencies C0=2n7lY where the oscillator is "frozen", i.e., < a >= 0. As we 
shall see this unusual behavior does not occur in our treatment because we always have a 
damped amplitude for all frequencies. 

The basic differences between our approach and Milburn' s approach (8) are: 
i) we do not make any specific statistical assumption: randomness in time evolution 

appears naturally as an interpretation of our results, i.e., of Eq.(3) and Eq.(8). In 

particular the Poisson-r distribution of Eq.(8) is not assumed but derived, 
ii) the choice Xi = X2 of ref. (8) implies a severe restriction on the statistical interpretation 

of the theory, 
hi) ref. (8) assumes that the number n of evolution transformations in a time t is a random 

variable and t is just a parameter as in the Schroedinger equation or in the Liouville 

equation. In our theory the waiting time t for n evolutions to occur is the random 

variable. 
The two view points appear similar but are basically different. In fact, they lead to 
different results, as one can see by comparing Eq.(14) and Eq.(7). The first can be 

derived from the second as follows. Using the integral identity 

-x 

ln(l + ix) = dX^— (l-e~ Ax ), and taking x = Lxi and X=t' /Xi, Eq.(7) can be written in 

* 
the integral form 



P~(t) = I 

X 2 V "l 




e- lLt '-l 



t , P(t) • (15) 

This equation is completely equivalent to Eq.(7). It can be easily shown that, if the term 
in square brackets is slowly varying on a time scale Xi, it can be taken out of the integral 
with t'= Xi. Because the term in the round brackets is normalized, one obtains the 
Milburn Master Eq.(14). Therefore, Eq.(14) can be obtained as an approximation of 
Eq.(7) or Eq.(15), if Xi=X2=Y _1 is small enough. In the case of an harmonic oscillator this 
implies the assumption co/y«1, which is inconsistent with the condition for freezing 
Cfl/Y=2n7C. 

CLASSICAL LIMIT AND THE GENERALIZED TAM 
MADELSTAM RELATION 



From the finite difference equation (10) and defining the mean value of an observable 
A in the usual way, A = Trp A , we obtain 

A(t)-A(t-T,) if 1 

x, - -»[*•"]• < 16 > 

Let us apply Eq.(16) to the one dimensional motion of a particle with Hamiltonian 

H = p 2 / (2m) + V(x) . Taking A=x and A=p x we have (x(t) - x(t - x 2 ) / x 1 = p x (t) and 

(p x (t) -p x (t-x 2 )/x 1 = F(x) , where F(x) = -dV/dx. This is the finite difference 
version of the Herenfest theorem. The usual form is regained in the continuos limit 
Tj = x 2 —>0, whereas the classical limit is F(x) ~ F(x) . Note that the two limits are 
independent. Therefore, taking only the classical limit one obtains finite difference 
Hamiltonian equations for x and p x . 

We now derive a generalized Tarn Mandelstam relation. Using the general uncertainty 



relation for A and H we can write a(A)o(H) > (1 / 2) [A, h] = ft|AA| / 2x 1 , where AA = 

A(t)-A(t-x 2 ) and we have used Eq.(16). The uncertainty relation can be written in 
the form x A o(H) > h/2;x A = o(A)/flAA \/x l J. This equation appears as a generalized 
TM inequality and reduces to the usual form in the limit Xi=X2— > 0. In fact, in this limit 
one obtains x A = o(A) / A . However, the uncertainty relation can be written as: 

IaaI x, 

L < — (17) 

o(A) x E V ; 

where Xe=^/2o(H) is the "intrinsic inner time" (14) of the system. Because inequality 

(17) is valid for any observable, we can conclude that if Xi<Xe, then AA <o(A) 

necessarily follows. Therefore, no appreciable variation occurs for any observable in the 
cronon time x 2 . We can say that if Xi<Xe one has a quasi continuous evolution of the 
system even using the discrete time description. On the other hand, Eq. (17) states that if 

%2 is such that AA > o(A), it follows that one must have Xi > Xe. Therefore, the ratio 

Xi/Xe rules the rate of change of the state of the system within the time interval x 2 between 
two evolutions. The choice Xi = Xe = h/2o(H) corresponds to the maximum possible 
value of Xi which guarantees the quasi continuous evolution which is commonly 
observed. This choice is clearly meaningless in case of a non Hamiltonian system. 



CALCULATION OF PHYSICAL QUANTITIES 

Equation (2) allows one to calculate all physical quantities by just performing a time 
integral of the usual expressions. Let us give some examples. The mean value 
A(t) = Trp(t)A can be easily obtained using by Eq.(2) as 

A(t) = jdt'P(t,t')<A> t , ; <A> t =Trp(t)A. (18) 

o 
Therefore, a constant of motion remains a constant of motion, whereas oscillating 

quantities, as e ltot will be damped with a rate constant y = (l/2x 2 ) ln(l + C0 2 x 2 ) similarly 

to p nm as given by Eq.(12). Furthermore, from Eq.(2), 

(x|p(t)|x') = jdfP(t,t')(x|p(t')|x') . (19) 



The same relation holds for matrix elements in any representation. However, Eq.(19) is of 
particular relevance. In fact, for x = x' it gives the position probability density, P(x,t) . In 

particular, if the initial state is a pure state, |\|/ ^\|/ |, then (x|p(t')|x) = |\|/(x,t')| . 

Furthermore, because the Wigner function is related to the Fourier transform of 

(x|p(t)|xy, it can be obtained using the same integral relation as in Eq.(19). It can be 

easily shown that the same relation is valid for the Glauber P-function. We now make 
some testable applications of our formalism. 

HARMONIC OSCILLATOR 

In the case of the harmonic oscillator the photon number is a constant of motion. On 
the contrary the amplitude <a> goes to zero. In fact, because < a > t =< a > e~ I<Bt , using 

Eq.(18) and Eq.(6) one obtains a(t) = (l/(l + iayc 1 ) t/T2 )< a > = < a > e _7t e" lVt , where 

Y = (l/2x 2 ) ln(l + G) 2 x 2 ) and v = (l/x 2 )arctgoyc 1 . Therefore, |a(t)| goes to zero as t 

goes to infinity as it must be for a phase-destroying process. The behavior of y and v is 
formally the same as that of y nm and v n>m described before. Therefore, the same discussion 
applies just abolishing the index n,m. This intrinsic linewidth y should be observable only 
if it dominates on all other linewidth, i.e., if Xi is large enough and % 2 is small enough. In 
the opposite case one obtains a lower limit for x 2 and an upper limit for Xi. 

FREE PARTICLE LOCALIZATION 

Let us consider a free particle prepared in a Schroedinger cat-like state i.e., in a 
coherent superposition state of two minimum uncertainty wave packets, centered at 
different positions, Xi, and x 2 . The wave function at time t is given by 

\|/(x,t) = (l/V2")(\|/ 1 (x,t) + \|/ 2 (x,t)) (20) 

where 



10 



V j (x,t)= '-—-^ 



-1 1 

°* J 



j=l,2. (21) 



(V2^o x ) 

Here, for simplicity, we have neglected the Schroedinger spread, by assuming o v to be 
small enough. This is certainly the case at the macroscopic limit, because 

o v =fil (2mo x ) , where m is the mass of the particle. From Eq.(20) we have: 

P(x,t) = |v ) /(x,t)| 2 =(l/2^ 1 (x)| 2 + (l/2^ 2 (x)| 2 +Int (22) 

where Int = |\|/ 1 (x)|\|/ 2 (x)|cos (Qt , with 

co = 16mE 2 Dx//z 3 , (23) 

where Xi = -x 2 = D/2 and E = (1/ 2)m < v 2 > . Therefore, in standard quantum mechanics 
the interference term is oscillating in time. In contrast, in our formalism the interference 
disappears because, using Eq.(18) and Eq.(6), the average of cos(u)t) is exponentially 

damped at a rate given by y = (l/2x 2 )ln(l + co 2 x 2 ), where (0 is given by Eq.(23). The 

interference disappears in a decoherence time to = Y~\ Note that if (QTi«l one can 
expand the logarithm to the first order, obtaining a decoherence rate proportional to (0 2 , 
i.e., to the square of the distance between the centers of the two packets, in agreement 
with previous treatments. The interaction with the environment and/or with the measuring 
apparatus appears quite naturally in our formalism via the two characteristic times 
without specifying any interaction model. Equation (23), together with the expression of 
y, clearly shows that at the macroscopic limit, when m and E are very large, or at the 
classical limit h — > the decoherence rate becomes extremely fast. 

INTRINSIC DECOHERENCE OF RABI OSCILLATIONS 

The same considerations can be applied to a system of two level atoms injected one at 
the time in a high Q resonant cavity and prepared in a Rydberg state so that the atomic 
and cavity decay time and the intrinsic decoherence are very long. This is the 
experimental situation in ref. (9). If the field is in a n-Fock state and the atoms are 
injected in the excited state, the atoms will oscillate between the upper and lower state so 

that the population difference d will oscillate asd = cos£2t, where Q. = gVn + 1 is the 
Rabi frequency and g is the one photon Rabi frequency accordingly with the Jaynes- 
Cummings atom-field interaction Hamiltonian model. In our formalism these oscillations 

will be damped with a rate constant, y = (l/2x 2 )ln(l + £2 2 x 2 ). Therefore, if yt = yL/v 

»1, the population difference will approach the steady state value d=0 so that the 
population of the upper level will approach the value 0.5. This behavior has been indeed 
observed in ref.(9) also for the vacuum state where dissipation losses are ineffective and 
in ref. (10) for a n-Fock state, in a different experimental situation. In particular, in 
ref.(10), it has been observed an increase of the damping rate y with n, in "qualitative" 
agreement with (n+1) ' 7 . Let us point out that our damping rate has a logarithmic 
dependence on (n+1), which, if Xi is small enough, goes like (n+1). Whether this damping 
is due to experimental problems, as suggested by the authors, or can be attributed to our 
intrinsic decoherence mechanism deserves further experimental and theoretical 
investigation. Obviously the same intrinsic damping can be predicted for the Rabi 
oscillations in a NMR configuration. 
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CANCELLATION OF EPR CORRELATION 

Let us assume two spin- 1/2 neutral particles, say 1 and 2, traveling in opposite 
directions with velocity v and prepared in the singlet entangled state |Vo) = 

(l/'v2j(|+ 1 ,— 2 ) — |— i,+ 2 )j with a constant magnetic field B in the z direction acting on the 
path of particle 1 for a length L. We have |\|/) = (l/V2J(e 1<Bot/2 |+,-)-e~ I<,)ot/2 |-,+)) where 

Cflb is the Larmor frequency. The associated density operator reads 

( 1 ^ 

p(t) = |i|/(t))(\|/(t)|= p D --e ifi, » t |+-)(-+|-h.c (24) 

V I ) 

where p D is the diagonal part. Note that for Cflt)L/v = 2n7l one has again the initial singlet 
state. Therefore, the presence of a e.m. field would become ineffective regarding EPR 
correlation under these conditions. However, in our formalism the oscillation of the off- 
diagonal terms disappears exponentially. Therefore, if yL/v »1 the EPR correlation 
should disappear even if 0>)L/v=2n7l. Any evidence of correlations smaller than the one 
expected in the conditions described above, would be a strong support of our theory. 

CONCLUSION 

We have generalized the Liouville equation without introducing new assumptions to 
quantum mechanics. On the contrary, we have reduced the basic axioms of quantum 
mechanics dropping the unitarity condition and maintaining only the semigroup property 
of the time evolution operator. The new equation describes intrinsic decoherence giving 
irreversible state reduction to the diagonal form in the energy basis. Several testable 
examples have been discussed. An intrinsic linewidth for a harmonic oscillator, space 
localization of a free particle prepared in a Schroedinger cat-like state and disappearance 
of EPR correlations under proper conditions have been described. In particular, damping 
of Rabi oscillations, in agreement with recent experimental results (9, 10) has been 
predicted. 
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